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ABSTRACT 

Black holes in core-collapse of massive stars are expected to surge in mass and 
angular momentum by hyper-accretion immediately following their formation. We here 
describe a general framework of extended emission in gravitational waves from non- 
axisymmetric accretion flows from fallback matter of the progenitor envelope. It shows 
(a) a maximum efficiency in conversion of accretion energy into gravitational waves 
at hyper-accretion rates exceeding a critical value set by the ratio of the quadrupole 
mass inhomogeneity and viscosity with (b) a peak characteristic strain amplitude at 
the frequency fb = Hb/vr, where Hb is the Keplerian angular velocity at which viscous 
torques equal angular momentum loss in gravitational radiation, with hchar oc at 
f < fb and hchar on at / > fb- Upcoming gravitational wave observations may 

probe this scaling by extracting broadband spectra using time-sliced matched filtering 
with chirp templates, recently developed for identifying turbulence in noisy time series. 


1. Introduction 

The birth of stellar mass black holes in core-collapse supernovae of relatively massive stars is 


(Bethe et a 

2nn.ii 

van Putten 

2nn4li 1 


extremal spin perhaps close to the Thorne limit (|Thorng|l974l L The black hole may subsequently 
experience spin-down by strong interactio ns with high density matter at the ISCO via an inner 
torus magnetosphere as (jvan Puttenlll999l L 


This general outlook on accreting stellar mass black holes is receiving increasing attention 
as a leading candidate of the central engines of long GRBs and low luminosity long GRBs and. 
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conceivably, the more extended class of broad-line hyper-ene r getic core collapse s upernovae (CC- 
SNe). Although an origin in millisecond magnetars ( Usov 1992 : Metzger et al. cannot be ruled 

out, rapidly rotating black holes offer multimessenger emissions through frame dragging induced 
interactions by an ample energy reservoir in spin energy, that can readily account for the most 
extreme GRB-supernovae (|van Putten et alJl2011a). By a diversity in spin, ro tating black holes 


further offer a unihcation of short and long GRBs ( van Putten fc OstrikenboOll ). 


The focus of the present discussion is hereby complementary to and different from that in 
existing models for gravitational waves from CG-SNe producing (proto-)neutron stars. The complex 
dynamical process of core-collapse and core-bounce in the first one or tw o seconds at birth of a 
neutron star opens a remarkably broad window to gravitational radiation (l0ttll2009ll of potentia l 
interest to upcoming observations of events at distances up to a few Mpc (jRover et ahl 120091 ) . 
The connection to GRB-supernovae and extremely energet ic supernovae, however, is not obvious 
( Burrows et al. 2007 : Dessart et al. 2008 : Mosta et al. 2015 1. 


Non-axisymmetric high-density accretion disks orbiting black holes are potentially prodigious 
emitters of gravitational radiation. A fraction of order unity of the gravitational binding energy 
of accreted matter released in gravitational waves would provide just the kind of energetic output 
sufficient for a detection within a distance on the order of H ~ 100 Mpc by advanced gravi- 


van Putten et al 


2911 , 

2Q14; 

van Putten 

2015 ) 


tational wave chirps from mass-inhomo geneities in accretion flo w, down to wave-instabilities about 


the ISGO in suspended accretion flows (jvan Putten 


2001 


2002) • The latter offers a window further 


to be powered by the spin energy of the black hole, giving a unique perspective on detecting black 
hole evolution. Ultimately, detection of gravitational wave emissions from a c ore-collapse super¬ 
nova is expect to rigorously identify the nature of the i nner engine of long GRBs (ICuttler &: Thorne 
2002), i-e., a black hole-disk system or magnetar (e.g. iMelatos fc PrivmakI (j2014l B. that presently 
remains elusive based on electromagnetic observations alone. 

Quite a few mechanisms can naturally give rise to non-axisymmetric accretion flows. First, the 
envelope of the progenitor is likely to be highly turbulent, leading to intermittent fallback of stellar 
material. Second, hydro dynamical instabilit ies give rise to unsteady accr etion with large temporal 


variations of the accretion flow parameters (jMacFadven &: Wooslevlll999r ). In 3D, it is conceivable 


that these instabilities excite non-axisymmetric azimuthal modes. A leading mechanism for gener- 


1965; 

Adams et al. 

1989 

Hadlev et al. 

20141. 1 

thin disks fe.e.. Gammie 

20Q1 

Rice et al. 

2005) 


becomes important. It c an even lead t o disk fraginentat i on if the cooliiig rate in the instability 
zone is sufficiently high (jGammiel 120011 : iRice et al.l l2005l : ISadowski et al.l l2015l ) . Although it re¬ 
mains to be demonstrated by direct numerical simulations, conceivably this mechanism applies also 
to collapsar disks at hyper-accretion rates relevant to long GRBs and some of the more energetic 
GG-SNe. Even without self-gravity, the accretion disk may also develop internal wave instabilities. 
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(Tagger et al. 

1990; 

Tagger & Pellat 

1999: 

2006; 

Lovelace & Romanova 

2014 

). Third 


20141 ') . Third, Papalo izou-Pringle wave instabilit i es may be excited 


at the surface of an inner torus about the ISCO (iPapaloizou Pringld Il984l : Ivan PuttenI l2002l i 
provided accretion is sufficiently suppre ssed (lBlaeall987li. e.g., by feedback from a rapid ly rotating 
Kerr black hole by magnetic coupling (jvan PuttenI 1 19991 : Ivan Putten fc LevinsonI 120031 ) . Fourth, 
accretion flows may be susceptible to transient shocks that m ay migrate over an extended ran ge of 
radii, ingoing or outgoing, with associated chirp-like behavior (jChakrabarti et al.ll2004l . l2008l l that 
conceivably break axisymmetry. Common to all these mechanisms is a gravitational wave output 
from mass-inhomogeneities in non-axisymmetic accretion flows. 


nection with black hole snin (van Putten 

2002 

• 1 

Putten & Levinson 20031. extended accretion 

flows ( 

Kobavashi & Meszaros 2003: 

Piro Sz Pfahl 

2007 

1, and recently also by numerical simulations 

(Korobkin et al. 

2011 

Kiuchi et. al 

2011: 

Mewes et al. 

20151. For self-gravitating disks, the Toomre 


criterion has been employed to argue that a collapsar disk should become gravitationally unstable 
at radii r/M > 300a^^fMi ( Piro fc Pfahl 2007 ). where a = 0.1a_i is the alpha viscosity pa¬ 
rameter and M = Mi Mqs~^ is the accretion rate. Several cooling processes have been examined, 
including electron-positron pair annihilation to neutrinos, URCA process, and photo-disintegration 
of '^He, with the latter fou nd to be the most effective one in the instability zone. It has been 
argued (jPiro fc Pfahlll2nn7l ) that cooling by photo-disintegration may lead to fragmentation chirps 
in gravitationally bound clumps, migrating inwards with frequencies 0.1 — 1 kHz depending on 
disk par ameters and fragmen tation mass. For LIGO, it suggests a sensitivity distance of about 
100 Mpc iPiro &: Pfahll (j2007l l. What fraction of the total mass accreted, if at all, will end up as 
gravitationally bound clumps, and how many clumps are formed are yet unresolved issues. 

Prodigious gravitational radiation at twic e the orbital frequency about the ISCO has been 
predicted in the suspended accretion scenario (Ivan Putten fc LevinsonI l2003l ) . In this model, the 
major fraction of the spin energy of a Kerr black hole is transferred to the surrounding torus via 
magnetic field lines anchored to it, and is emitted as gravitational waves, MeV neutrinos and winds. 
A minor fraction is released along open magnetic field lines in the polar region in the form of a 
ultra-relativistic jet. The torus is subject to a self-regulated instability, with a quadrupole mass 
inhomogeneity of the order of 10% at saturation level. It is worth noting that in this scenario, the 
energy source of the gravitational waves emitted by the torus is the rotational energy of the black 
hole, rather than the gravitational binding energy of the surrounding matter, as in the gravitational 
radiation from accretion flows. The suspended accretion model requires a st rong suppression of 
mass loading of magnetic field lines that connect the torus and the black hole (j Globus fc Levinson 
20131 ). possibly supported by surface gravity in the inner face of the torus by it s inherently super - 
Keplerian motion. It could be that this process may be inherently intermittent (jvan Puttenlll999l l. 
Observational support for the associated loss of angular momentum of the black hole, following 
a near-extremal state formed in a prior epoch of hyper-accretion, is found in normalized BATSE 
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light curves (|van Puttenll2012l . 120151 ') 


In light of the various distinct physical processes by which non-axisymmetric accretion flows 
may originate in core-collapse events, we here consider a general framework for their gravitational 
wave emissions in hyper-accretion flows following black hole formation. Our primary objective is 
to identify some canonical scaling behavior in broadband emissions from radially extended non- 
axisymmetric features in accretion flows. Our focus is on relatively long-lasting features potentially 
lasting tens of seconds in accretion flows derived from fall back matter from a progenitor stellar 
envelope, posterior to a relatively brief but complex transition in the formation and settling down 
of the inner accretion disks immediately foll owing black hole formatio n, recently highlighted in 
detailed axisymmetric numerical simulations ([Sekiguchi &: Shibatall201lh . 


Specihcally, we set out to identify the frequency at maximal gravitational wave luminosity 
following a transition between viscous and gravitational wave-driven accretion flows. These general 
results may be conhrmed in future large scale numerical simulations of accretion flows, that take 
into account gravitational wave emission and their back reaction on the accretion flow. 


2. GWs from non-axisymmetric accretion flow 


We consider a mathematical discretization of radially extended non-axisymmetric features in 
accretion disks by annular rings, each rotating at the local Keplerian angular velocity 


f2(r) 


1 

M VT J ’ 


( 1 ) 


around a black hole of mass M. We henceforth use geometrical units, G = c = 1, whereby mass 
accretion rate is measured in units of 

c3 

hio = -^ = 4 X 10^8 g s-\ (2) 

and luminosity in units of Lq = mgc^. The underlying non-axisymmetric features themselves may 
be continuous, as in spiral waves, or discrete, in case of lumps. In what follows, we shall develop 
a general frame work for the associated broadband gravitational wave emission based on our ring 
discretization, that applies irrespective of the detailed properties these features. The condition for 
our ring discretization to apply is rather mild, merely that it identihes mass-inhomogeneities 6m 
therein concentrated over a limited azimuthal extend, i.e., 6ipj2Tr « 1 in each ring of sufficiently 
small width 6r (Fig. 1). To be precise, consider a spiral wave pattern in the orbital plane (x = 
r{ip)cosif,y = r{if)siinp) in polar coordinates (r, y?) satisfying \dr{ip)/dip\ > k > 0 for all ip. 
(Here, k = 0 is excluded, as dr{ip)/dip = 0 corresponds to a circle.) Then 6(p < St/k, and hence 
bipj^-K << 1 for br sufficiently small. 

We next consider a Lagrangian description, in which the accretion disk is partitioned in rings. 
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A ring located at a radius r has a mass 


Am{r) = a{r)M 


(3) 


and a radial width /(r). It is naively anticipated that Z(r) is on the order of the vertical scale height 
H(r) of the disk. For a thin ring of radius r rotating at a frequency Q in the (x, y) plane, the 
non-vanishing components of the quadrupole moment tensor, 

= f cfxp{t,x)x^x^, (4) 

Jv 

where the integration is over the ring’s volume, are given, to a good approximation, by 

jxx _ ^ cos 2Qt), 

jyy = iAmr^(l — ^ cos 2nt), (5) 

I^y = — sin 

in terms of the dimensionless parameter ^ that quantihes the quadrupole mass inhomogeneity, 
whereby ^ = 0 for an axi-symmetric ring, ^ = 0.5 for a ring having a density p{r, p) = po cos^ p, 
and .^ = 1 for two identical point masses in circular orbit around their center of mass. 

For simplicity, the dynamics of the ring and its gravitational wave emission are computed in 
the Newtonian limit. This is a reasonable approximation far enough from the ISCO. However, as 
shown below, the gravitational wave power is maximized near the ISCO, where absorption and 
scattering of gravitational waves and other GR effects are at play. Thus, our analysis should be 
treated as a first approximation. 


Orbital angular momentum is lost both in viscous torques and gravitational radiation, leading 
to an inward radial drift of the ring. For > 0 the gravitational wave power of the ring, ALcwi^)^ 
is assumed to be dominated by Z = m = 2 mode described by the quadrupole formula. 


dEcw _ 

dt 5 dP’ dt^ 


where PP = i.e.. 


QO 

ALcwP) ^ 

5 




(7) 


With our assumption on m = 2 dominance, we focus on rings that are luminous in gravitational 
waves and mechanisms of non-axisymmetric instabilities for which low order modes tend to be 
unstable first. Furthermore, the gravitational wave output of a collection of rings that are each 
m = 2 unstable will be the sum of ([7j), each at its own gravitational wave frequency defined by 
twice the Keplerian frequency (fT|). 

An illustrative example is a spiral wave pattern in Fig. [TJ Overlaying our discretization of 
the plane of the disk in annular rings, this pattern introduces over-dense regions in each ring. 
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Each ring hereby assumes a m = 2 quadrupole moment inherited from the underlying spiral wave, 
here represented by a pair of mass-inhomogeneities 6m in each ring. Due to rotation, each 6m 
emits quadrupole gravitational wave emission at twice the local Keplerian frequency. Subject to 
accretion, the tidal wave gradually tightens. As the 6m gradually migrate inwards, these quadruple 
emissions broaden in frequency, each taking their own window in the gravitational wave spectrum. 
By Parseval’s theorem, the total luminosity from the spiral wave is the sum of the luminosities in 
each ring. 


It may be appreciated that our focus on quadrupole gravitational wave emissions is a leading 
order approach, that ignores the possibility of emissions from higher order mass moments. A 
formalism for calculating the full spe ctrum of gravitational wave emission fro m multiple mass 
moments in an annular ring is given in lBromberg. Levinson fc van PuttenI (j200fl). 


For gravitational wave driven accretion, the loss of orbital angular momentum in gravitational 
waves causes a radial displacement dr of the ring on a time 


dtcw ~ 


d{AE) 

ALgw 


5 



^ ^cT(r) ^dr, 


( 8 ) 


where AE = —Am/2r is the energy of th e ring in the N ewtonian limit. Viscous torques acting on 
the ring cause a radial drift dr on a time (jPringldll98lh 


dtu 


2rdr 

3u 


3a 



/ r \ 1/2 

\m) 


dr, 


(9) 


where the alpha prescription is adopted to express the kinematic viscosity in terms of the vertical 
scale height H{r) and the a parameter as u = aCgH = aVlH^/y/2, using c* = QH/y/2 between the 
sound speed and angular velocity D. The two times are equal at r = r;, given implicitly by 


rl'^Hl I23y/2 il 

M5/2(J^ 15 a ’ 

here the subscript b denotes values at r^. 


( 10 ) 


For viscous driven accretion, i.e., at radii where angular momentum losses are dominated by 
viscous torques, the surface de nsity of the disk is given approximately by S(r) = m/(37rz^) for 


Pringld (jl98ll ). The width of the ring can be identified with the scale height of 


accretion rate m 

the disk, hence we set l{r) = H(r). Then, the mass of the ring can be approximated as 


cj(r) = 27rrZ(r)S(r) 


Vs 

3a \H 


r \3/2 

mJ 


Combining Equations (fTOll and (fTTI) yields 


( 11 ) 


El 

M 


45 o? 


m. 


El 

n 


2 


( 12 ) 
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At r > r;, viscous stresses contribute the dominant torque, whereas at risco < r < r u the torque is 


dominated by gravitational radiatio n. Semi-analytic models of neutrino cooled disks iPopham et al 


( 1999lhlChen k, Beloborodov! ( 20071 1 indicate that H/r changes from about 0.1 near the ISCO to 0.2 
at r = lOM to 0.45 at r = lOOM, with relatively weak dependence on a and M for the conditions 
relevant to our analysis. We shall therefore adopt for illustration H(r) = rjr, rj = 0.1r7_i, in the 
inner disk region. Denoting o-i = a/O.l, we obtain a rough estimate for the transition radius: 


a. 

M 


Tifi ~ 






(13) 


We thus arrive at m and (fT^ constraining rf,. 

In the outer region r > r^ the mass profile is given by Equation (HU). Combining with Equation 
([7|) yields 

AT / \ f 

ALcwir > n) = ~ ^ ( 1 ^) 

45 V ^ / 

In the inner region risco < r < ri, the mass profile of the rings is determined from the relation 
adr/l = rndtcw- Then, by employing Equation ([8|) one has 

< rt) = , (15) 

and ALgw{^ < fb) = rhr]M/2r, upon substituting Equation (fT^ into Q. 

The total luminosity, Lqw^ is the sum over the contributions of all the rings, each radiating 
at their own frequency as indicated above. Approximating the sum by an integral we can write: 
Lew ~ ALGw{'f)dr/I. In the parameter regime of inefficient emission, 

n < '^iscoi (16) 


the mass profile given by Equation (HU holds everywhere and hence, using Equation HU , we 
further have 


Lgw r°“‘ 256 i'^M'^dr . 128 < ^ 

Mc^ Jrisco '^5 a‘^H^ 45 


(17) 


i.e., 

^^15^i^ ( M \ 

Mc^ yi Mq s ^ J 


(18) 


for rout » risco- Here < >= r ^^‘^{r)dr, and z = risco/M. The parameter regime 


^ '^isco 

is different. In this regime Equations (fT5|) and ([7|) yield 

Lgw ^ J_ 

Mc2 “ 2z 


(19) 


( 20 ) 


















for ri) » Tisco- Efficient emission occurs at accretion rates for which Equation (jl9p is satisfied. 
With aforementioned definition of z, Equation (jl3p gives the corresponding criterion 

M > 0.157?^ 1 >-i a^zMQS-^ (21) 

whenever a mass inhomogeneity assumed to originate at larger radii, survives all the way to the 
ISCO. If it dissipates at some radius r^iss > fisco, then risco should be replaced by rdiss-, thus raising 
the minimum accretion rate on the right hand side of (|21l) . 


Eor a matched filtering detection method, the relevant quantity is the amplitude that takes 
into account the square root of the associated number of wave periods. In the frequency doma in, 
the corresponding quantity is the characteristic strain amplitudel (e.g., ICuttler &: Thornd 120021 '): 


hcharif^ 


ttD 


AE 


A/ 


( 22 ) 


Consider now coherent emission episode of duration r. In the accretion model outlined above, r 
corresponds to the period during which a steady disk pattern, excited by some instability, survives, 
satisfying 


tew < T <T, (23) 

where T is the total accretion time. The energy radiated over time r from an annular region located 

between r and rout » r, assuming ri, > r > rtsco, is 

front /M\ ( r \ 

AErad = r J ALGw{r')dr' ~ mr f — j il - —j . (24) 

We remark Eqn. (j24p is independent of that governs the gravitational wave luminosity Eqn. ©• 
As such, ^ governs the time scale over which AE^ad is realized. A small (large) ^ has a low (high) 
luminosity in GWs, whereby AE^ad takes a long (short) time to come to fruition. The net result, 
however, is always the same AErad- 


The energy AErad is emitted over the frequency interval Af = /(r) — f{rout) = ^{r)/Tr — 
Q{rout)/'^ — 7r“^M“^(M/r)^/^, for the assumed Keplerian rotation. Thus, the characteristic strain 

amplitude is Khar if > h) = ^\ I , i-e., 


A/ 


dchar{,f ^ fb) 



M ^ ^ 



(25) 


denoting fb = fKb). In the outer region r > rb we employ Eq. (fT^ to obtain AErad — tALgw oc 
{r/M)~^, emitted over the frequency interval Af ~ fir)- Thus, hcharif < fb) oc Careful 

calculations yield: 




if < fb) 

ifisco > f > fb) 


f^charif) — ^ ' 


1/6 


(26) 
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with K scaled to distances D 2 = (D/lOOMc) and masses Mi = M IOMq satisfying 


= V8 




rriT 


O.IMq 


1/2 


h 


lOOOHz 


-1/6 


(27) 


Here fisco = firisco) is the frequency in quadrupole GW emission of matter at the isco. Eq. (I26]l 
applies also in the regime fb > fisco, whereby hchar = K{f/fb)^/^, f < fisco- 


Fig. [2] displays a plot of Eq. (f26]l for different values of fb and mr. It shows the characteristic 
strain amplitude peaks at the frequency fb, given by 


fb 





M 


1 Mq s 


-1 


-3/2 


Hz, 


(28) 


where Eq. m has been used. The condition > rj^co implies fb < fisco, where for a IOMq black 
hole 430 Hz < fisco < 3000 Hz, depending on the spin parameter a/M of the black hole, with 
fisco = 1600 Hz for a/M = 0.95 as an example. 


Fig. [3] shows the frequency range relevant to gravitational wave emission from hyper-accreting 
black holes; Fig. [^exhibits a plot of fb versus M, elucidating the dependence on the spin parameter. 


3. Conclusion and Outlook 


Fig. [ 2 ] indicates a window of extended gravitational wave emission, that may arise from long- 
lived quadrupole mass inhomogeneities in hyper-accretion flows. As a result of accretion, these 
inhomogeneities will produce chirps, as opposed to constant frequency emission such as may be 
produced by a pulsar. In deriving our scaling relation (j28|) for the peak luminosity in gravitational 
waves, as illustrated in Fig. 2, we display a radiatively efficient and inefficient parameter regime 
(fT6l) and, respectively, (fT9]) . depending on the viscosity and accretion rate in the accretion flow. 


Broadband emission in chirps can be searched f or broadband spectra ex tracte d from noisy time 


van Putten et al. 

(2011 

))■ In 

van Putten et al. 


(j2f)14l )l. some 8.6 million templates were used to densely cover the intended range in frequency 
and frequency time rate-of-change. In the application to gravitational wave searches, exhaustive 
searches of this kind require massively parallel computing and a relatively accurate determination 
of the True Time Of Onset (TTOO) of the supernova, which may be estimated by extrapolation of 
i ts optical light c urve, that may be obtained from dedicate supernova surveys of the Local Universe 
( Heo et al. 2015 1. 


A LIGO-Virgo or KAGRA detection of the proposed broadband emissions promises entirely 
novel constraints on the viscosity and accretion rates in core-collapse supernovae, not currently 
accessible by existing electromagnetic observations. Considerable gain in sensitivity can be attained 
by coherent network analysis of the four detector sites of these observatories (e.g. IWen fc Schutz 
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(|2005l l: iHavama et alJ (j201a)). Detection promises unique insight into the different phases of birth 
and evolution of stellar mass black holes. Probes of broadband extended gravitational emission 
such as illustrate d by spiral waves in Fig. [1] may, for instance, identify a critical accretion rate for 
hyper-accretion (jGlobus fc: Levinsonll2014l b the end of which leaves a near-extremal black hole as 
the initial condition to p ronounced IS CO w aves during prompt GRB emission in subsequent phase 
of black hole spin down (jvan Puttenll2ni5l ). 
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Fig. 1.— (a) A spiral density wave pattern in a disk, (b) Identification of over-dense regions 
(thick black) in an annular region 4.5 < r < 5 (thin black circles) that are of finite angular extend 
(5y?/27r << 1. (a) Leading order approximation of the over-dense regions in (b) by local mass- 
inhomogeneities 6 m, here 2 x 32 in number following a grid with 32 annular rings, (d) By their 
individual and distinct Keplerian angular frequencies, the quadrupole emission spectra of each 6 m 
in (c) are non-overlapping. Expressed in terms of the modulus squared of the associated Fourier 
coefficients Ck, accretion broadens the emission of each, shown in (d) for two rates of accretion 
corresponding to different radial migrations 6 r 2 < 6 ri < 0. 
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Fig. 2 .— Characteristic strain amplitude hcharif) of broadband quadrupole emission in accretion 
flows around a stellar mass black holes of M = IOMq at a fiducial distance oi D = 100 Mpc. 
The vertical distance to the dimensionless strain noise hn = VfSh of LIGO S5 (2007) represents 
the theoretical limit of signal-to-noise ratio in a matched filtering detection. The different lines 
correspond to different values of fh and Ma = mr, as indicated. Thick lines correspond to the 
maximally efficient parameter regime fb < fisco in Eq- thin lines correspond to ft > fisco- 
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Fig. 3.— Schematic overview of quadrupole gravitational radiation from hyper-accreting black 
holes. Wave patterns in accretio n flow may produce broadband emission (circles), whereas 
fragmentation (Piro fc Pfahll 2007 ) or non-axisymmetric waves in a torus about the ISCO 
( van Putten &: LevinsonI 2003) (red curve) may produce chirps. Relatively high frequency emis¬ 
sion may derive from Quasi-Normal Mode (QNM) ringing in the black hole event horizon (black 
curve). 
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Fig. 4.— A plot of ff, versus accretion rate M, for the canonical choice of parameters adopted 
in Eq. (j28p . The horizonthal lines delineate the ISCO frequence for two choices of the spin 
parameter, as indicated. The region enclosed above the fh line and below the fisco line (diamonds) 
indicates the regime of maximum efficiency of gravitational wave emission. The circles indicate the 
permitted frequency range of less efficient emission, though still relevant for detection by advanced 
gravitational wave detectors, as indicated by the thin lines in Fig. [2l 





